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We study the all-optical switching behavior of one-dimensional metal-dielectric photonic crystals 
due to the nonlinearity of the free metal electrons. A polychromatic pump-probe setup is used to 
determine the wavelength and pump intensity dependence of the ultrafast transmission suppression 
as well as the dynamics of the process on a subpicosecond timescale. We find ultrafast (sub- 
picosecond) as well as a slow (millisecond) behavior. We present a model of the ultrafast dynamics 
and nonlinear response which can fit the measured data well and allows us to separate the thermal 
and the electronic response of the system. 

PACS numbers: 



I. INTRODUCTION 

In t he la te 1980's, the field of photonic crystals 
emerged^El. Since then It has rapidly evolved and 
grown more interesting to a wide community Photonic 
crystals show remarkable properties due to their band 
structure, which lead to the realization of a number of ap- 
plications, such as bending of light around 120° corners^, 
extremely slow light in a medium^ and negative refrac- 
tion while keeping the index of refraction positive^. 

Another promising field of application is the combi- 
nation of photonic and electronic devices, or even the 
replacement of electronics by photonics or plasmonics 7 . 
Such all-optical circuits have several advantages over con- 
ventional electronics, such as reduced size, high repeti- 
tion rates, and enhanced speed of operatiorP. 

A key component for such circuits is an all-optical 
switch. One-dimensional metal-dielectric photonic crys- 
tals (lDMDPCs) might fulfill this purpose, as small 
pump-induced variations in the dielectric properties can 
lead to strong transmission changes. In this article we 
study experimentally the suitability of lDMDPCs for ul- 
trafast switching in the subpicosecond regime, determine 
the detailed wavelength and intensity dependences and 
give a theoretical model, which agrees well with our ex- 
perimental data. 

A 1DMDPC is a structure which exhibits a periodic 
change of the index of refraction in one dimension. Al- 
ternating layers of metal and dielectric are deposited with 
a subwavelength thickness. This gives rise to a pho- 
tonic bandstructure with passbands and stopbands in 
the transmission spectrum. A large amount of metal can 
thus be accumulated in these structures while keeping 
them transparent. This unique feature allows to exploit 
the nonlinear properties of the constituting materials 10 . 
Various properties and possible applications have been 
proposed based on numerical calculations and experi- 
ments, for e xample op tical limiting, switching and selec- 
tive shieldin g* 9 ! 11 ! 12 ! 13 *. However, the time constant of the 
nonlinear respose of such structures was not measured. 
While previous work made use of the strong nonlinearity 
of the bound metal electrons, we explore the possibility 



of using the response of the free electrons. That approach 
has the advantage that no material resonance is involved 
and thus the whole device could be tailored over a much 
larger spectral range. 



II. EXPERIMENTAL TECHNIQUES 

We use a polychromatic pump-probe setup to investi- 
gate the 1DMDPC. The laser system consists of a Ti:Sa 
oscillator, which is pumped with 10 W at 532 nm. This 
gives rise to 1.5 W average power and a pulse length 
of 200 fs at 825 nm. The repetition rate of the mode- 
locked laser system is 76 MHz. We split the beam into 
a pump and a probe arm. The pump arm length is tun- 
able mechanically, which allows high temporal resolution 
of our setup, limited only by the pulse length. We mod- 
ulate the pump beam with an acousto-optical modulator 
(AOM) at 100 kHz for lock-in detection. In a second 
beampath, a tapered fibei^ is used to generate a white 
light continuum from which the probe pulse is shaped in 
the time as well as in the frequency domain. This al- 
lows us to choose the probe pulse central wavelength and 
width while still keeping the pulse duration at 600 fs. 
The pulses are overlapped in a polarizing beamsplitter 
cube and focused collinearly on the sample. Behind the 
sample the pump beam is filtered out spectrally, and a 
balanced receiver detects the transmitted probe power. 

The samples are prepared using an electron beam evap- 
oration technique. Alternating layers of Ag and MgF 2 are 
evaporated on a quartz substrate at a pressure of 10 -7 
mbar. The metal is evaporated with a rate of 2QA/s to 
avoid clustering and to achieve smooth films, since oth- 
erwise silv er ten ds to cluster when applied at such small 
thicknesse d 15 ! 16 ! The smoothness of the film, on the other 
hand, is critical for consistent measurements. For reasons 
that become obvious below, we refer to a layer series of 
metal dielectric metal as a cavity. Two samples were 
used in this work: a single cavity sample (30 nm Ag and 
225 nm MgF2) and a double cavity sample (20 nm Ag 
and 215 nm MgF2). The first has a narrow transmission 
resonance at 720 nm, the latter at 705 and 810 nm. 
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FIG. 1: center: Measured linear transmission spectrum of the 
'double cavity' 1DMDPC. Passbands with multiple transmis- 
sion peaks occur, marked A and B. The simulated spectrum 
was calculated using the transfer matrix method. The red dot 
indicates the fixed pump wavelength used in this paper, left 
and right: The electric field distributions at the wavelengths 
marked A and B superimposed on the structure (Ag: yellow, 
MgF2: purple). The field forms standing waves, with a sym- 
metric (B) and anti-symmetric (A) mode and nodes only at 
the metal layers. The transmission peaks at smaller wave- 
lengths correspond to higher order modes with an additional 
node in the center of each MgF2 layer. 



III. LINEAR TRANSMISSION PROPERTIES 

In a first step, it is instructive to look at the lin- 
ear transmission spectra of lDMDPCs and their band- 
structure. Figure [l] shows a measured transmission spec- 
trum of the single cavi ty sam ple. Calculations using the 
transfer matrix method^^ reproduce the features very 
well. Each passband (around 700-850 nm and around 
300-450 nm) is approximately 150 nm wide and con- 
sists of two peaks. The small third peak in the pass- 
band at 310 nm is due to the interband transition of sil- 
ver. To unders tand these features, we consider the Bragg 
conditio n 1 19 ! 20 !. The metal layers form Bragg planes, 
which are separated by dielectric layers of thickness d. 
This leads to standing waves inside the structure, if the 
condition n ■ d = k ■ Xq/2 is fulfilled, where n is the refrac- 
tive index, Ao is the wavelength in vacuum and k is an 
integer. However, the Bragg model is only a first approx- 
imation, as according to it the long wavelength passband 
should occur at 675 nm and have only one peak. 

We can also regard our system as coupled cavities. 
Each layer series of metal-dielectric-metal forms a cav- 
ity. Consecutive cavities are coupled by the electric field 
leaking through the joining metal layer. A series of N 
cavities can be described by a series of N coupled os- 
cillators which exhibit N eigenmodes. Figure [T] shows 
the electric field distribution for the symmetric and anti- 
symmetric eigenmode of the lowest-order passband, cal- 
culated with a FDTD method. The second-order modes 
(around 300-450 nm) have an additional node of the field 
in the center of each cavity. 




FIG. 2: Open aperture z-scan measurement of the double 
cavity sample at a pump wavelength of 825 nm (see red dot in 
Fig.flJ. The sample was measured with a continuous wave and 
a pulsed laser beam with the same average power (20 mW), 
but peak powers differing by about five orders of magnitude. 
The inset shows the calculated differential transmission due 
to thermal effects. The assumed temperature modulation is 
AT = 50 K and 500 K starting from room temperature. 



The separation between two eigenfrequencies, i.e., 
two transmission peaks, is given through the coupling 
strength of the oscillators, i.e., the thickness of the silver 
layer. A thinner metal layer corresponds to a stronger 
coupling of the cavity system, since the evanescent light 
field decays exponentially inside the silver. The trans- 
mission is lower for the symmetric fundamental mode 
(marked B), as the electric field undergoes a ir phase 
jump at the silver layer, which leads to a higher field 
inside the metal and therefore to a greater loss. 



IV. THERMAL EFFECTS 

In contrast to experiments at a low laser repetition 
ratelSl, our repetition rate of 76 MHz leads to a larger 
power absorbed in the sample at even moderate pulse 
energies. The temperature rise can reach several hun- 
dreds of Kelvins locally, and thermal effects have to be 
considered. Any change in the optical properties due to 
temperature will lead to a cavity detuning and thus to 
a change in transmission which could mislead the inter- 
pretation of experiments that do not contain temporal 
information on t he picosecond timescale, such as z-scan 
experiments^^!] A z-scarP^ is a single beam experi- 
ment, where the sample is moved through the focus along 
the optical axis (z-axis) of the beam. The transmission 
is measured depending on the z-position of the sample, 
i.e., as function of the laser intensity, as the focus size 
varies. Such measurements are used to determine the 
nonlinear absorption coefficient and the nonlinear refrac- 
tive index. Our laser system allows us to switch between 
mode-locked and continuous wave mode, changing thus 
the peak power by about five orders of magnitude without 
changing the average impinging power. Figure [2] depicts 
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the results of an open aperture z-scan measurement of the 
double cavity sample at a pump wavelength of 825 nm 
and a focal length of 50 mm. In our case, the change in 
transmission is not due to effects that are nonlinear in 
the optical field strength, but is only related to the de- 
posited energy in the system and thus the heatload. The 
thermal effect has the expected slow rise and decay times 
(230 ms and 1570 ms, respectively). 

A thermal response of the 1DMDPC can have several 
sources: thermal expansion due to the thermal expan- 
sion coefficient 5, change of the index of refraction due 
to the thermo-optic coefficient T, change of the index of 
refraction due to the nonlinear refractive index ri2, and 
the change of the plasma frequency of the metal in an 
expanded layer. The last two points are negligible, since 
the layer expansion and thermo-optic effects are orders 
of magnitude larger for the materials used in our experi- 
ment. The inset in Fig.[2]shows the calculated differential 
transmission for temperature variations AT of 50 K and 
500 K starting from room temperature using the transfer 
matrix method. The parameters^ for this calculation 
are S Ag = 1.89 • 10~ 5 K~\ 8 Mg F 2 = 1 ■ 10~ 5 K~ l , and 
TivigFa = 2 • 10~ 6 K~ x . The calculated effect reaches the 
same order as the measured transmission change of the 
z-scan. 



V. TIME-RESOLVED ULTRAFAST 
EXPERIMENTS 

To avoid the ambiguity of z-scan measurements with 
respect to thermal effects, we now turn to time-resolved 
pump-probe experiments on the picosecond timescale. 
They will allow us not only to separate the influence of 
heat accumulation on a millisecond timescale from ultra- 
fast effects, but also to gain the time constants of the 
ultrafast response of the 1DMDPC. 



A. Ultrafast response of metals 

The ultrafast response of the 1DMDPC is governed by 
the response of the metal layers, wh ich can be separated 
into four phases that overlap in tim d 24 ! 25 ^: 

- Pump photons excite some free electrons into higher 
levels far above the Fermi level, which leads to a non- 
Fermi-Dirac distribution, since most of the electrons are 
still 'cold'. Our pulses of 600 fs length can not resolve 
this short phased 

- In a few tens to hundreds of femtoseconds, the elec- 
tron gas thcrmalizcs via electron-electron scattering to 
a Fermi-Dirac distribution. Electron gas temperatures 
of several hundreds to a few thousand Kelvins can be 
reached. As our experiment in the near infrared probes 
the free electrons, this phase leads to the highest system 
response. 

- The electron gas transfers its energy to the lattice 
via electron-phonon coupling. Due to the large heat ca- 



pacity of the lattice, a single laser pulse increases the 
temperature by only a few ten Kelvins. 

- The coupled system cools down to the same tempera- 
ture as outside the laser spot. The accumulated absorbed 
energy together with the slow heat conduction away from 
the laser spot leads to the temperature increase of some 
hundreds of Kelvins, as discussed in the last section. 

As soon as a thermal distribution of the excited elec- 
trons is reached, the coupled electron-lattice system can 
be described by the two-temperature model (TTM^PU. 
The temperatures of the electron gas and the lattice (T e 
and Ti) are given by two coupled differential equations: 

8T - - 

C e ^- = -g(T e -T l )+X7K\7T e + S (1) 

C l -^-=g(T e -T l ) . (2) 

C e = aT e is the specific heat capacity of the electron gas, 
k is the heat diffusion coefficient, S is the source term, 
g describes the electron-phonon coupling, and Ci is the 
specific heat capacity of the lattice. The temperatures 
and the source term are functions of space and time. To 
simplify the calculations, we avoid a full spatial treatment 
of the system by neglecting the spatial dependence and 
by replacing the exact heat diffusion term with an empiric 
diffusion term including a fit parameter d, which couples 
the system to a temperature bath at room temperature 
Troom- This accounts for the heat transport out of the 
system, which is hereby assumed to be governed by the 
electron gas. The equations are then reduced to 

dT e _ g (T e - T{) S 
dt ~ aT e d[le lmom)+ aT e [6) 

cm = g(T e -Tj) 

dt d 1 } 

The parameters^- 4 for the simulations are a = 
66J/(m 3 K 2 ), g = 3 ■ l0 16 W/(m 3 K), and Q = 2.415 • 
10 6 J/(m 3 K). We chose Ci to be constant, since we mea- 
sure at a temperature well above the Debye-temperature 
of silver (ss 200 K) . The source term is assumed to be 
a Gaussian pulse (in time) which is absorbed in a spot 
of 10 fj,m diameter. For reasons of simplicity, we assume 
the whole pulse energy to be absorbed in these numerical 
calculations. 

For the calculation of the transient transmission sup- 
pression we incorporate the TTM into the transfer matrix 
and reference the work of Bigot et oZffl The electron gas 
and lattice temperature are calculated using the TTM. 
To account for the change of the optical properties, the 
dielectric function e(ui) of the metal layer is modified, 
starting from Johnson and Christy values^, with a Drude 
term for the free electrons that depends on the rise of the 
electron gas temperature AT e : 

e(ui, AT e ) =e J+c (w) - e D rudc(w, AT e = 0) 

+ eDrudo(w, AT e ) . (5) 



4 



The Drude model is made temperature dependent by let- 
ting the damping constant 7 depend on AT e . This is mo- 
tivated as with increasing electron gas temperature the 
Fermi distribution smears out and more and more possi- 
bilities for electron-electron scattering occur. We setP^l 



7 (AT e ) = 7o + /3-AT e 



(6) 



where (3 is a fit parameter in our calculations and 70 the 
room temperature value of the damping constant. The 
modified Drude model reads thus 



£Drude(^, AT e ) — £ c 



^7, 



lu{u + i 1 (AT e )) 
where ui p is the plasma frequency. 



(7) 



B. Temporal dynamics of the nonlinear optical 
response 

Our pump-probe setup (Fig. [3^,) allows us to investi- 
gate not only the amplitude but also the temporal dy- 
namics of the nonlinear optical effects in lDMDPCs. A 
pump pulse is absorbed, modifies the dielectric proper- 
ties of the metal layers, and a probe pulse interrogates 
the transmission of the photonic crystal after a variable 
delay. Figure [3J) shows a series of such measurements on 
the single cavity sample with varying intensity. The dif- 
ferential transmission is not zero for negative delays (i.e., 
probe arrives before pump) because of the thermal effect 
discussed above. This offset decreases with increasing 
AOM modulation frequency, as less heat is accumulated 
in each AOM cycle. The probe transmission drops at 
the arrival of the pump pulse and recovers on a subpi- 
cosecond time scale. It does not recover fully to the start 
value, since also the last pump pulse deposits additional 
heat in the focal area. 

To extract quantitative data from these traces and sim- 
ilar ones from the double cavity sample, we fitted a phc- 
nomcnological model function /(r) to the data: 



f{r)=A + 



Be 



-t/t 



C(l 



-t/tq\ 



9(t) 



(8) 



A describes the amplitude of the thermal offset at neg- 
ative delays. B is the amplitude of the response caused 
by the hot electron gas that decays exponentially with a 
time constant tq. The thermal effect of the pump pulse is 
taken into account by a level C that is exponentially ap- 
proached with the same time constant tq. The last two 
processes start as soon as the pump pulse arrives, i.e., 
the Heaviside step function 0(t) switches from (r < 0) 
to I (t > 0). To take into account the limited temporal 
resolution of our setup, we convolute the model function 
/(t) with a measured cross-correlation function between 
pump and probe pulse before fitting it to the data. A 
fit to a temporal trace at high pump power is shown in 
figure [Sk 
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FIG. 3: (a) Detection part of the experimental setup, 
(b) Pump-probe measurements with varying pump intensity 
(given in MW/cm 2 as numbers in the graph) on the single 
cavity sample. The absolute zero delay was set to the differ- 
ential transmission peak minimum. The sample was pumped 
at 825 nm and probed at 720 nm. (c) Example of the model 
function (Eq. [8| fitted to one data set from (b). 



C. Intensity dependence of the nonlinear optical 
response 

The amplitude B and the decay time To of the non- 
linear optical response were recovered by fitting the phe- 
nomenological model function /(r) to a set of temporal 
traces with different pump intensities. Figure|4]shows the 
data for the single and the double cavity sample. The ab- 
solute value of the nonlinear optical response \B\ seems to 
grow linear with pump intensity for both samples, always 
maintaining its negative sign. The decay time To is not 
constant but also increases with pump intensity. Both ef- 
fects can be described by a TTM that modifies the Drude 
damping in the metal's dielectric properties as described 
in section |VX Only two fit parameters are used for both 
parts of Fig. |4j firstly the relation between the increase 
in electron gas temperature AT e and Drude damping 7, 
described by /? in equation ([6]), and secondly the empiri- 
cal constant d to describe the heat diffusion by electrons 
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FIG. 4: The TTM combined with a temperature-dependent 
Drude model describes with one parameter set pump-power 
dependence of both (a) the peak depth (amplitude B) as well 
as (b) the recovery time constant to of the ultrafast transmis- 
sion suppression. The latter is weakly dependent on the pump 
intensity due to the temperature dependence of the electron 
heat capacity. 



in equation d3). The fit leads to = 6 • 1CT 8 eV/K and 
d = 3.7 • 10 1/s. The deviations at small pump in- 
tensities for the single cavity sample in figure |4j) might 
occur due to a lesser sample quality and to slightly differ- 
ent probe locations on the sample, which lead to locally 
changed heat diffusion and therefore to changed dynam- 
ics. 

The temporal response in the lDMDPCs is on the 
order of 1 ps. The effect of rising time constants with 
increasing pump intensity can be understood, since the 
heat capacity of the electron gas is proportional to its 
temperature (C e = aT e ). The lattice heat capacity 
is orders of magnitude larger, but higher pump inten- 
sities lead to a smaller difference between electron gas 
and lattice heat capacity. This leads to a slower equili- 
bration (longer response time) and to a slightly higher 
end temperature. When the electron gas heat capacity 
is assumed to be constant, for example at its value at 
T e = 400 K (C e = 2.64 • 10 4 J/(m 3 K)), then the time 
constant is independent of the pump intensity (see blue 
line in figure |4Jd). This confirms the temperature depen- 
dent electron heat capacity as origin of the pump inten- 
sity dependence of the decay time. The general trend of 
increasing response times at higher pump intensities was 
also found in single layers of metaPS. The values of the 
time constants were on the order of 1-2 ps. 
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FIG. 5: Spectral dependence of the transient transmission 
compared to the linear transmission spectrum and the model 
for (a) the single cavity sample and (b) the double cavity sam- 
ple. The pump power (1.3 GW/cm 2 ) and pump wavelength 
(825 nm) was kept constant. The TTM simulations give a 
good agreement. 

To compare the magnitude of the nonlinear optical 
response of the lDMDPCs with other materials, it is 
instructive to calculate the phcnomenological nonlinear 
absorption coefficient a e ff- It describes the intensity de- 
pendence of the absorption coefficient a 

a(I) = a Q + a e gl. (9) 

As the absorption coefficient a is connected to the trans- 
mission T via the sample thickness L: T — Toe~ aL , the 
nonlinear absorption coefficient a e ff is just a measure for 
the slope of a linear fit to the data in figure HJi. We ob- 
tain values of 5.2 • 10~ 10 cm/W and 7.4 • W~^cm/W for 
the single and double cavity sample, respectively, using 
as sample thickness L the sum of all metal and dielectric 
layers in the multilayer structure. 

D. Wavelength dependence of the nonlinear optical 
response 

Figure [5^ shows the wavelength dependence of the ul- 
trafast transmission suppression in the single cavity sam- 
ple. For this measurement, pump power and pump wave- 
length were kept constant (1.3 GW/cm 2 at 825 nm). The 
shape of the differential transmission data follows the lin- 
ear transmission spectrum, which means that the spec- 
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tral response of an all-optical switch could be tuned by 
adjusting the thickness of the metal and dielectric layers. 

The double cavity sample shows a similar behavior 
(Fig. [5j)). The TTM simulations predict a stronger re- 
sponse at the 810 nm transmission peak compared to the 
peak at 705 nm. The same trend is found in the exper- 
iment. However, the probe spectral range is limited by 
the pump filter stopband starting at 780 nm. The mea- 
surement shows a steeper rise than the simulations at 
775 nm. This might be due to the roughness and inho- 
mogeneity of the silver films. The position of the 810 nm 
peak is dependent on the silver layer thickness, and the 
position of the 705 nm peak is relatively insensitive to 
thickness variations. The white light transmission spec- 
trum is measured with a large spot (4 mm in diameter), 
which leads averages over the roughness. The pump- 
probe measurements, on the other hand, were performed 
with a spot size of roughly 10 fim. If a region with a 
slightly larger silver thickness was measured, the result 
would be a shift of the 810 nm peak into the blue, which 
would explain the rise at smaller wavelengths. The TTM 
parameters for these simulations are f3 = 1.1 • 10~ 8 eV/K 
and (3 = 3.5 • 10~ 8 eV/K for the single and the double 
cavity sample, respectively. 



infrared pump pulse train at a Megahertz repetition rate: 
on a timescale of hundreds of milliseconds thermal effects 
lead to a large transmission change of 30-40%. Simu- 
lations based on the transfer matrix were presented to 
explain the measurements. On a picosecond timescale, 
the ultrafast response of the metal's free electrons dom- 
inates. Its magnitude is in the order of 10~ 5 in a sin- 
gle cavity system and 10 -3 in a double cavity system. 
The two-temperature model and a temperature depen- 
dent Drude model explain both the pump intensity and 
the probe wavelength dependence of the effect. Using 
the nonlinear response of the free metal electrons makes 
it possible to tune the spectral response of the metal- 
dielectric photonic crystal by its structure, i.e., without 
changing the material. However, that flexibility is accom- 
panied by a red uced nonlinearity compared to material 
resonance d 13 * 21 [ 
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